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ABSTRACT 

We present an approximate solution to the minimally coupled Einstein-Dirac equations. 
We interpret the solution as describing a massive fermion coexisting with its own gravita- 
tional field. The solution is axisymmetric but is time dependent. The metric approaches 
that of a flat spacetime at the spatial infinity. We have calculated a variety of conserved 
quantities in the system. 



1 Introduction 



Explicit solutions in General Relativity (coupled to matter) and in supergravity theories 
play a significant role in the study of quantum gravity. Several important progresses (for 
example, [H El O S] ) have been achieved with the help of some concrete solutions. Due to 
the weakness of the gravitational interaction, most of the phenomenologically interesting 
solutions describe physics at the astronomical scale, such as black holes, p-branes, black 
rings [5] and so on. 

From the theoretical point of view, however, it is also interesting to look for exact 
solutions at the microscopic scale. Quantum field theory (QFT) in flat spacetime tells 
us that if the zero-point energy gravitates, then extreme fine tuning is needed to render 
the cosmological constant to its currently observed value. There have been a lot of effort 
towards solving the problem, but without too much success [6]. The zero-point energy is a 
byproduct of the quantization of fields in flat spacetime. When it comes to curved spacetime, 
the situation is much more complicated (see, e.g. 0). As implied by the Unruh effect [2], 
even the notion of particle is no longer fundamental, but depends on the observer. It is 
our hope that an exact solution to gravity coupled to quantum matter may help us better 
understand the nature of QFT in curved spacetime, and offer some clue to the persisting 
problem of the cosmological constant. When a spacetime is flat at the spatial infinity, 
the usual notion of particle is still useful for an observer far away from the distribution of 
matter. One can imagine a scenario where the observer sees nothing but a single particle 
at the center. It is then interesting to ask how the wave function of the particle (such as a 
neutrino or electron) behaves under its own gravity. To find this, we will also need an exact 
solution to classical gravity coupled with quantum matter. More recently, the development 
in non- Fermi liquid and holographic superconductors [H [9l [1^ has also made it interesting 
to study gravitational solutions with back reaction from spinor fields [11]. In this work, 
we will present an approximate solution to the Einstein-Dirac system. It describes how a 
fermion exists, under its own gravity, in a flat spacetime background. 

Compared with the long list of literature on other systems, there have been relatively 
fewer work on solutions to the Einstein-Dirac system (for examples, see [12j-|26j). Most 
of the effort came after Brill and Wheeler's 1957 paper [12]. An early review of works 
focused on neutrinos can be found in |13j . Some of the solutions describe ghost spinors (see, 
e.g. [H]). Solutions in dimensions other than four can be found in [15\ 116] . Particle-like 
solutions have been previously studied in [25] . 

Apart from the complexity of the calculation, the reason discouraging people from deal- 
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ing with the spinor field might be its anti-commutative nature upon quantization. Here we 
will treat the spinor field as a quantum mechanical wave function, but not as a quantized 
field operator. The reason is the following. For a bosonic field, the wave-function treatment 
can be justified in the low energy limit, when there is no particle production. On the other 
hand, a spinor field is often believed to be intrinsically quantum, and one should always 
quantize it first before using it. However, the problem with gravity is that there is NO 
known ways to quantize the spinor field that is completely satisfactory. There is always 
a negative energy associated with the un-quantized spinor field. In the usual QFT treat- 
ment, this negativeness is shifted to hide in the zero-point energy, rendering it physically 
irrelevant [27] . But in theories with gravity, the problem reappears in the form of a too 
large contribution to the cosmological constant. In this sense, the cosmological constant 
problem is essentially the negative energy problem. So when there is gravity, the usual 
way of quantizing the spinor field is no longer a good fix of all the problems. Our hope is 
that, an exact particle-like solution to the Einstein-Dirac system (but with an un-quantized 
spinor field) may offer some hint on how to correctly quantize the spinor field with gravity. 

Due to the complexity of the coupled Einstein-Dirac equations, we start by looking for 
an approximate solution to the system. We expand the solution in terms of the radial 
coordinate, and our solution is approximately valid at the spatial infinity. At the present 
stage, it is very difficult to discuss the stability of the solutionis] In fact, it is not even clear 
if our approximate solution comes from a well behaved full solution or not. On the other 
hand, it is already non-trivial to find an ansatz (such as ()12p and ()13p ) which allows for 
an approximate solution to exist So as a first step towards finding an exact particle-like 
solution to the Einstein-Dirac system, we will be content with the approximate solution for 
now, but will try to solve all the remaining problems in future works. 

I will present the solution in next section. Then I will discuss some of the (approxi- 
mately) conserved quantities. A short summary is at the end. 



2 The Solution 



The action of the Einstein-Dirac system is given by 




R-2A i 



(1) 



thank the referee for raising up this point. 
^For example, the counterpart of (|12p and (|13p allow for no solution when there is an arbitrary cosmo- 
logical constant. 
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where is Newton's constant in n spacetime dimensions, A is the cosmological constant, 
and fi is the mass of the spinor field. The signature of the metric is mostly positive. It is 



often convenient to use the Planck mass M. 



P ~ SttG 



(n) 



in place of Newton's constant. We 



will let 167rG(„) = 1 from now on. The notations related to the spinor field are 



a ; 



Note 7*^ in -f/J is defined in the vielbein basis. The equations from ([TJ are 



(2) 



R 



(3) 



We will work in four spacetime dimensions {n = 4) from now on and will take A = 0. As a 
matter of choice, the gamma matrices in the vielbein basis are taken to be 



7° = i 



7l'2.3 = i 



a' 



3,1,2 



-a 



3,1,2 



(5) 



where I2 is the two-dimensional unit matrix and o"^'^''^ are the usual Pauli matrices. The 
corresponding charge conjugating operator is C = 7^. 

We are interested in an axisymmetric spacetime. The ansatz for the metric is given by 



ds^ = -f^{dt + f,d(t>f + fldx^ + fdO^ + {fpdcj) - f.dtf 



(6) 



where \ft\ > \fz\ in general. We expect x to be the asymptotic radial coordinate, 9 the 
latitudinal angle, t the time and (j) the azimuthal angle. It is convenient to write the 
vierbeins as 



= ftidt + hdcP) , e^ = hdx, = fydO , = fpdcj) - f.dt . 
One can always set = by using the following local Lorentz transformation 

/ 1 fz/ft\ 

1 n 



(7) 



1 - ifz/ft)' 







\fz/ft 



(8) 



1 / 



So ([6]) and ([7]) are equivalent to the following, 



.0 _ 



e- = Mdt + /, 
We take the spinor field to be 



fxdx . 



fyde, 



fpdc/). 



(9) 
(10) 



(11) 



1p2a + ilp2b 
tp3a + ilpSb 

where all the functions ipia and ipib (through out the paper, the index i = 1, • • • ,4) are real. 

Given ([9]) and pT]) . it is difficult to solve ([3]) and @ directly. In this work, we will focus 
on a spacetime that is flat at the spatial infinity (x — t- +oo). Our strategy is to expand the 
functions in terms of the radial coordinates x and obtain an approximate solution to the 
equations when x — t- +oo. We will assume that all the unknown functions depend on x,^ 
and t only. Our ansatz for the functions are 



ft-- 

fx 

fp 



l + ^ + ^ + 0(^), /. 

1 



£l 

X 



£2 
2 



+ ^ + o(- 



X 



x-^ 



1 



IH ^^ + C'(^), /y = x + yoH V^ + 0[-^) 



sm ( 



x + po + - + ^ + 0(A) 



Oil , , 1 N , , ^j2 , 1 ^ 



I,--- ,4: 



(12) 

(13) 



where xi,2,... , ii,2, - , ^1,2, - , ?/o,i,2,- , Po,i,2, - , Oji,j2, - and 6ii,j2, - are polynomials of 
trigonometric functions of x, and i. The expansion in (I12p guarantees that ([9]) approaches 
the metric of a flat spacetime as x — ?• +00. 

Given (fT2]) and (fT3]) . the leading contribution to ([3]) is of the order 0{^). The equations 
are given by 



- {dx - dt)aii = 
f^hi + {dx + dt)a2i = 
l^ibu + {dx + dt)asi = 
I^b2i - {dx - <9t)a4i = 



/iflsi + {dx - dt)bii = 0. 
^041 - {dx + dt)b2i = 0, 
mi - {dx + <9t)63i = 0. 
/^a2i + {dx - dt)b4i = 0. 



(14) 



These equations can be completely solved by 

flii = ajia(0) cos(t(;t) cos(/cx) + aji{,(0) cos(t(;t) sin(/cx) 
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+aiic{0) sm{wt) cos{kx) + andiO) sin(wt) sin(A:x) , (15) 

where i = 1, • • • ,4 and w'^ = k'^ + /u^. feji's are solved in terms of a^i's in an obvious way. 
The leading contribution to (jH) is of the order 0{x). The equations are given by 

(dl - d^)yo = (dl - df)po = . (16) 

They do not depend on the spinor field. The solutions are 

yo = yoy{6) + yoi{0,x-t) + yo2{9,x + t), 

PO = POy{0)+Poi(.0,X-t)+po2(.0,X + t). (17) 

It is hard to see how can functions of wt it kx and xzizt coexist with each other in the case 
w ^ k. So we will assume yoi = 2/02 = Poi = P02 = when fJ- ^ (In fact, our solution 
exists only when ^ 7^ 0). Note the leading order equations (jl4p and (I16p are independent 
of each other. 

To the order O(^) of ([4]), one has equations that determine zi,ti,xi,yi and pi, in 
terms of an and bn. The assumption about the functions zi,ti,xi,yi and pi impose extra 
constraints on the structure of an and bn. For example, two of the equations are 

\aub21dx In (^] + a2ibudx In 



dxdyti 



dydtxi 



4yr^ 



yb2iJ ^ Kb 
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4x/T-r ^ 



-031641 In (^) - a^ibsidx In (^) > (18) 
aiib2idt In f + a2ibiidt In 



,2 



-O3i64i0t In — - 0416319* In 

V 641 / V 631 



(19) 

Using (lisp , one can find that the right hand side of (llSp does not depend on x, while that 
of (I19p does not depend on t. To be consistent with the assumption that ti and xi only 
depend on j; or t in the form of trigonometric functions, the right hand sides of both (llSp 
and (I19|) must vanish. This condition constrains both ti,a;i and anjbn. As a special case 
that satisfies this condition, we let 

611 ~ 021 , 621 ~ ail , 631 ~ 041 , 641 ~ 031 , (20) 



up to some constant factors. This will reduce (jl5p to a much simpler form. 

The equations from higher orders of the expansion are more complicated and we will 
only summarize some of the main features here. There are eight equations at each order 
0{l/x^) of ([3]). These equations are just enough to determine a^n and bin in terms of 
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Ui^n-i, bi,n-i, Xn-1, Un-i, Pn-1, tn-1, Zn-1 and lower Order functions. On the other hand, 
the equations at the order ©(l/x") of (jH determine Xn-i, Vn-i, Pn-i, tn-i and in 

terms of at^n-l, h^n-l, Xn-2, yn-2, Pn~2, t 

n—2i Zn—2 and lower order functions. What's 
more, there are always two equations that involve Zn,tn,Xn,ymPn, o-in and bin at each order 
0{l/x^) of dH), while the rest only involve lower order functions. These equations can always 
be solved together with those from the order of Because of the constraints 

from equations like (jlSp and ()19p . it is not guaranteed that one can find a self-consistent 
solution to all orders in i. We have done the calculation up to the order 0{l/x'^) for both 
([3]) and dH), and still we do not see any true obstacle (other than the tediousness) to push 
the calculation to higher orders. This is an encouraging sign that (I12p and (I13p might be 
the right ansatz to give us a consistent solution. 

At the moment, we will be content with the solution approximate up to the order 
0{l/x'^). Even at this stage, the full result is already very unwieldy and contains many 
free parameters and functions. Presumably these free functions and parameters should be 
determined by equations from higher orders of the expansion. But for the purpose of giving 
an accessible example, we have (not so rigourously) chosen the parameters and functions 
in such a way that the solution has a simpler structure. One cannot expect to get the 
correct approximation to some exaction solution (if it exists) in this way. But the solution 
so obtained is still approximately valid in its own right, in the sense that it solves ([3]) and 
dH) up to the order ©(l/x^). 

For the metric, we find that 



ft 



X 



y 



X + 



1 - 



1 - 




IGk'^wfi'^x'^ ' 
iV2 \quj cos(2A;x) sin^ O + k sm{2kx) sm{2wt - 2C) sin(26') 



+ 



Sk'^fi'^x'^ 




(21) 



For the spinor field, we find 



cos(w;t + A:x — C) 



iV2 



(4A; — w)w cos cos{wt + kx — Q 



a 




Sk'^lix 



7 



+ 



{2k + w) cot 9 sm{wt — kx — () 



lb 



1p2a 



1p2b 



3a 



'4a 



'4b 



4:k{w + k)x 



•^^/w + k - 



cosi^wt — kx — Q) 
+ 



N'^{4:k + w)w cos 9 cos{'wt — kx — Q) 



Sk'^fix 

{2k — w) cot 9 sm{wt + kx — (^) 
4:k{w — k)x 



+ 0{l/x^ 



smt/ r {Ak — w)w cos 9 cos{wt — kx — C) 
cos{wt — kx — 



•^^/w — k - 



+ 



Sk'^fix 

{2k + w) cot 9 sm{wt + kx — () 



4:k{'w + k)x 



+ 0{l/x^ 



•:y/'w + k 



cos{wt + kx — 
4 



N'^{Ak + w)w cos 9 cos{wt + kx — () 



Sk'^fix 

{2k — w) cot 9 sm{wt — kx — Q) 



4k{w — k)x 



+ 0{l/x^ 



smt 



z\Jw — k - 



sm.{wt — kx — Q 
+ 



N'^{Ak — w)w cos 6 sm{wt — kx — C,) 



Sk'^jix 

{2k + w) cot 9 cos{wt + kx — 



4k{w + k)x 



ArVski^ 



z^/w + k - 



sm{wt + kx — Q) 
+ 



N'^{4:k + w)w cos 9 sin(u;t + /cx — C) 



Sk'^fix 

{2k — w) cot 9 cos{'wt — kx — Q) 



4:k{w — k)x 



+ 0{l/x^), 



smt 



z\Jw — k 



sm{wt + kx — Q 
+ 



N'^{4:k — w)w cos 6 sm{wt + kx — Q 



Sk'^fix 

{2k + w) cot 9 cos{wt — kx — Q) 



4k{w + k)x 



+ 0{l/x^), 



A^VskT^ 



Zy/w + k - 



sm{wt — kx — Q) 
+ 



N'^{4k + w)w cos 9 sin(u;t — kx — Q) 



Sk'^fix 

{2k — w) cot 9 cos{'wt + kx — C) 



+ 0{l/x^), (22) 



4k{w — k)x 

where w'^ = k"^ + , N is a. dimensionless normalization constant and C, is an arbitrary 
phase. For the results given above, it is possible to set C = by a shift in t. The function 
cot 9 appearing in ()22p indicates that the solution is divergent at sin = 0. This could be 
a big problem. But such divergence could also be an artifact of the expansion in (jl2p and 
psp. For example, something like (a and h are constants) 



b + x sin 9 

is obviously regular at sm9 = 0, but it diverges at sm9 = if firstly expanded around 
x — )• +00. For (j22p . it is possible to get rid of the divergence in a similar fashion. 

We have found ([2T]) and ([22]) by assuming that a^n and bin are polynomials of {cos(/cx). 
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sm{kx)} and {cos{wt), sin(wt)} up to the (2n — l)'th power, and Zn, tn, Xn, Vn and pn are 
polynomials of {cos{kx), sin(A:x)} and {cos(wt), sin(^yt)} up to the 2n'th power. The depen- 
dence on 9 appears as coefficient functions of the polynomials. This assignment is inspired 
by the structure of (jlSp . (jl7p and that of the equations from different orders of ([3]) and (j3]). 
After solving equations up to the order 0{l/x'^), we find that many coefficient functions 
are still undetermined. What's more, we have also found a lot integration constants in the 
process. Although these functions and constants are arbitrary at the order 0(1/2;^), one 
should expect many further constraints to arise from higher order equations in the expan- 
sion. So presumably these functions and constants should be determined by pushing the 
calculation to higher orders, which is currently a daunting task. In reaching (|2ip and (I22p . 
we have thrown away most of the functions that do not depend on the strength of the spinor 
field (i.e. A^) in an obvious way. These functions may be important when one wants to go 
to higher orders. As a result, (12ip and (|22p may not be the correct approximation to an 
exact solution of ([3]) and ([4|). On the other hand, (I2ip and ()22p is still a valid approximate 
solution, good to the order 0{l/x'^). 

3 Conserved Quantities 

From ([3]), one can derive a conserved current of the spinor field, 

j^^ = ^^f^^^ =^ D^J^' = 0. (23) 

If there is a time- like killing vector ^, then one can write down the probability density of 
the spinor field as 

P = J-^- (24) 

For the energy density of the spinor field, there are two possibilities. One is by using the 
energy-momentum tensor, 

Pi = eeT^u , = R^u - ^g^u . (25) 

The other is by using the energy- momentum four- vector, 

P2 = P-L P^ = -li^D^i; + ^D^iPip. (26) 

It can be checked that V^P^ = 0. We will see that pi and p2 are quantitatively vastly 
different. 
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From ()2ip . one can find an approximate time-like Killing vector ^ = 5t at the 

spatial infinity x — )■ +00, 

d?yi\ 



( 



dgdtyi 



dgdtyi 

X 



2 SVC? 9 dty 2 



\ X 



(27) 



where /U, 1/ = {x,6,(j),i}, and yi,y2 can be read off (fT2]) and ([2T]l . Now the probability 
density of the spinor field can be calculated as 



M?; AN'^w sin ( 



2 



+ o(^). 



(28) 



Here we have restored the Planck mass to make the dimension of the density manifestly 
correct. If we interpret the solution as describing a particle with its center of mass located 
at X = 0, and if we suppose that the chance of finding the particle far away from x = is 
g, then 

2tt^M^N'^wL 



dVp 



where 



dV 



dx / d9 







2n 



sin ( 



(29) 



(30) 



The integral (j29p is divergent over the whole space, so we have introduced a cutoff L to 
regularize the divergence. Now if the full wave function is finite near the center, then we 
must have q = \. This is because the distribution probability is then dominated by the 
divergent integral in ([29l) . We will keep q explicit to cover the possibility that the wave 
function may actually diverge at x = 0. In this case, a significant portion of the probability 
may be distributed near the center. 

Let's now turn to the energy of the spinor field. The energy densities (f25]l and (f26|) are 
found to be 



92 



/\/f2 N^w COS 9 



2/i + SiV^/c sin(2'u;t - 2C) sin 6 



2 2A;/x2x3 



+ 0(l/x^ 



(31) 



Here we have again restored the Planck mass to make the dimensions manifestly correct. 
It is obvious that the two energy densities are very different. While the nature of p\ is not 
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very clear, we will take p2 to be the true energy density of the spinor field at places far 
away from the center. In fact, it is easy to see that 

P2 = pii + 0{l/x^). (32) 

The contribution of p2 to the total energy is 

E = j dVp2 = Qp . (33) 

In the case q = 1, p2 dominates the contribution to the spinor energy. We see that the 
total amount is exactly p. This is as expected: the solution describes a particle without 
an apparent kinetic energy, and so the total energy is nothing but the mass of the particle. 
However, there is a puzzle. In the solution the wave function fluctuates in both space and 
time, and the frequency w is apparently larger than the mass. We must have a non- vanishing 
wave number {k ^ 0) for the solution to exist. But it is still unclear how k is related to the 
particle mass p. One may need to know the full solution to answer this question. 
On the geometry side, the Ricci scalar of (|2ip is 



^^2iV^.^sin(2..t-2C)sm.^ (34) 
px-^ 

It is obvious that the curvature can be both positive and negative, and oscillates with a 

frequency of 2w. Since the spacetime is flat at the spatial infinity, one can use the Komar 

formula to calculate the energy stored in the geometry [28l [29], 

^ r M^Ti'^kN'^L 

M = Ml = ^ sin(2wt - 2C) 

^ Jx=L fJ- 

= --^sm{2wt-2C). (35) 
2 w 

This energy fluctuates in time as well and it also goes negative for half of the time. We do 
not have a good explanation to this result at the moment. For an answer, one may have to 
better understand how gravity is coupled to quantum matter. 

Similar to ()35p , one can try to calculate the angular velocity of the geometry, 

J = / *d{d^) = ^-^ cos(2fcL) « 0. (36) 

^ Jx=L "J/^ 
It will be interesting to compare this result with the angular momentum of the spinor field. 

By analogy with the result in a flat spacetime, we look at the quantity 

^'' = ^eV^W"V', with \e^upa\ = V\9\- (37) 

We find that all the components vanishes as in p6p . If this result is also true for the full 
solution, then it means that the particle described by the solution does not have a fixed 
axis of spin. 
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4 Summary 



We have presented an approximate solution to the Einstein-Dirac system. It solves the cou- 
pled Einstein-Dirac equations up to the order 0{l/x'^), with x being the radial coordinate. 
The solution can be interpreted as describing a single Dirac fermion coexisting with its own 
gravitational field. The metric approaches that of a flat spacetime at the spatial inflnity. If 
one assumes that the full wave function is everywhere regular in the whole space, then the 
total energy in the spinor held is just the mass of the particle. The energy in the geometry 
fluctuates in time, and it is negative for half of the time. For the solution to exist, we also 
need a non-vanishing wave number k in the radial direction. The value of the wave number 
is undetermined and we still know very little about its significance. 

A natural generalization of the present work is to include a cosmological constant in the 
spacetime background. We will leave this to future works. 
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